** Replication of Fear and Loathing by Brooklyn Walker and Donald P. Haider-Markel, published in Public Opinion Quarterly

* Use Fear and Loathing Replication Data

* CN = Christian nationalism
* BA = born again
* Ideo = ideology
* SR = symbolic racism
* WhDisc = perception of discrimination against Whites
* CDisc = perception of discrimination against Christians


** Table 1: Effect of Treatment

reg Worried i.Treatment 
reg Angry i.Treatment
reg Afraid i.Treatment
reg Disgusted i.Treatment

** Table 2: Indirect Effects of Treatment on Christian Nationalism

* Race treatment
sem (Disgusted -> CN, ) (Angry -> CN, ) (Afraid -> CN, ) (Worried -> CN, ) (PartyID -> CN) (Ideo -> CN) (SR -> CN, ) (SR -> Disgusted, ) (SR -> Angry, ) (SR -> Afraid, ) (SR -> Worried, ) (BA -> CN, ) (BA -> Disgusted, ) (BA -> Angry, ) (BA -> Afraid, ) (BA -> Worried, ) (Attend -> CN, ) (Attend -> Disgusted, ) (Attend -> Angry, ) (Attend -> Afraid, ) (Attend -> Worried, ) (Edu -> CN, ) (Income -> CN, )  (Male -> CN, ) (Age -> CN, ) (TreatRaceControl -> CN, ) (TreatRaceControl -> Disgusted, ) (TreatRaceControl -> Angry, ) (TreatRaceControl -> Afraid, ) (TreatRaceControl -> Worried, ), vce(bootstrap, rep(500) dots(1))  cov( e.Disgusted*e.Angry e.Disgusted*e.Afraid e.Disgusted*e.Worried e.Angry*e.Afraid e.Angry*e.Worried e.Afraid*e.Worried) nocapslatent

estat eqgof
estat gof, stats(all)

nlcom _b[Afraid:TreatRaceControl]*_b[CN:Afraid]
nlcom _b[Worried:TreatRaceControl]*_b[CN:Worried]
nlcom _b[Angry:TreatRaceControl]*_b[CN:Angry]
nlcom _b[Disgusted:TreatRaceControl]*_b[CN:Disgusted]

nlcom [Angry]_b[TreatRaceControl]*[CN]_b[Angry] + [Worried]_b[TreatRaceControl]*[CN]_b[Worried] + [Disgusted]_b[TreatRaceControl]*[CN]_b[Disgusted] + [Afraid]_b[TreatRaceControl]*[CN]_b[Afraid]

* Religion treatment

sem (Disgusted -> CN, ) (Angry -> CN, ) (Afraid -> CN, ) (Worried -> CN, ) (PartyID -> CN, )  (Ideo -> CN, )  (SR -> CN, ) (SR -> Disgusted, ) (SR -> Angry, ) (SR -> Afraid, ) (SR -> Worried, ) (BA -> CN, ) (BA -> Disgusted, ) (BA -> Angry, ) (BA -> Afraid, ) (BA -> Worried, ) (Attend -> CN, ) (Attend -> Disgusted, ) (Attend -> Angry, ) (Attend -> Afraid, ) (Attend -> Worried, ) (Edu -> CN, ) (Income -> CN, )  (Male -> CN, ) (Age -> CN, ) (TreatRelControl -> CN, ) (TreatRelControl -> Disgusted, ) (TreatRelControl -> Angry, ) (TreatRelControl -> Afraid, ) (TreatRelControl -> Worried, ), vce(bootstrap, rep(500) dots(1)) cov( e.Disgusted*e.Angry e.Disgusted*e.Afraid e.Disgusted*e.Worried e.Angry*e.Afraid e.Angry*e.Worried e.Afraid*e.Worried) nocapslatent

estat gof, stats(all)
estat eqgof

nlcom _b[Worried:TreatRelControl]*_b[CN:Worried]
nlcom _b[Disgusted:TreatRelControl]*_b[CN:Disgusted]
nlcom _b[Afraid:TreatRelControl]*_b[CN:Afraid]
nlcom _b[Angry:TreatRelControl]*_b[CN:Angry]

nlcom [Angry]_b[TreatRelControl]*[CN]_b[Angry] + [Worried]_b[TreatRelControl]*[CN]_b[Worried] + [Disgusted]_b[TreatRelControl]*[CN]_b[Disgusted] + [Afraid]_b[TreatRelControl]*[CN]_b[Afraid]
display   .597 / (.597+.467 )

** White discrimination DV

* Race treatment

sem (Disgusted -> WhDisc, ) (Angry -> WhDisc, ) (Afraid -> WhDisc, ) (Worried -> WhDisc, ) (PartyID -> WhDisc, )  (Ideo -> WhDisc, )  (SR -> WhDisc, ) (SR -> Disgusted, ) (SR -> Angry, ) (SR -> Afraid, ) (SR -> Worried, ) (BA -> WhDisc, ) (BA -> Disgusted, ) (BA -> Angry, ) (BA -> Afraid, ) (BA -> Worried, ) (Attend -> WhDisc, ) (Attend -> Disgusted, ) (Attend -> Angry, ) (Attend -> Afraid, ) (Attend -> Worried, ) (Edu -> WhDisc, ) (Income -> WhDisc, )  (Male -> WhDisc, ) (Age -> WhDisc, )  (TreatRaceControl -> WhDisc, ) (TreatRaceControl -> Disgusted, ) (TreatRaceControl -> Angry, ) (TreatRaceControl -> Afraid, ) (TreatRaceControl -> Worried, ), vce(bootstrap, rep (500) dots(1))  cov( e.Disgusted*e.Angry e.Disgusted*e.Afraid e.Disgusted*e.Worried e.Angry*e.Afraid e.Angry*e.Worried e.Afraid*e.Worried) nocapslatent

estat eqgof
estat gof, stats(all)

nlcom _b[Worried:TreatRaceControl]*_b[WhDisc:Worried]
nlcom _b[Disgusted:TreatRaceControl]*_b[WhDisc:Disgusted]
nlcom _b[Afraid:TreatRaceControl]*_b[WhDisc:Afraid]
nlcom _b[Angry:TreatRaceControl]*_b[WhDisc:Angry]                                                                                             
nlcom [Angry]_b[TreatRaceControl]*[WhDisc]_b[Angry] + [Worried]_b[TreatRaceControl]*[WhDisc]_b[Worried] + [Disgusted]_b[TreatRaceControl]*[WhDisc]_b[Disgusted] + [Afraid]_b[TreatRaceControl]*[WhDisc]_b[Afraid]
display .063 / (.063+.196)

* Religion treatment

sem (Disgusted -> WhDisc, ) (Angry -> WhDisc, ) (Afraid -> WhDisc, ) (Worried -> WhDisc, ) (PartyID -> WhDisc, )  (Ideo -> WhDisc, )  (SR -> WhDisc, ) (SR -> Disgusted, ) (SR -> Angry, ) (SR -> Afraid, ) (SR -> Worried, ) (BA -> WhDisc, ) (BA -> Disgusted, ) (BA -> Angry, ) (BA -> Afraid, ) (BA -> Worried, ) (Attend -> WhDisc, ) (Attend -> Disgusted, ) (Attend -> Angry, ) (Attend -> Afraid, ) (Attend -> Worried, ) (Edu -> WhDisc, ) (Income -> WhDisc, )  (Male -> WhDisc, ) (Age -> WhDisc, ) (TreatRelControl -> WhDisc, ) (TreatRelControl -> Disgusted, ) (TreatRelControl -> Angry, ) (TreatRelControl -> Afraid, ) (TreatRelControl -> Worried, ), vce(bootstrap, rep(100) dots(1)) cov( e.Disgusted*e.Angry e.Disgusted*e.Afraid e.Disgusted*e.Worried e.Angry*e.Afraid e.Angry*e.Worried e.Afraid*e.Worried) nocapslatent

estat eqgof
estat gof, stats(all)

nlcom _b[Worried:TreatRelControl]*_b[WhDisc:Worried]
nlcom _b[Disgusted:TreatRelControl]*_b[WhDisc:Disgusted]
nlcom _b[Afraid:TreatRelControl]*_b[WhDisc:Afraid]
nlcom _b[Angry:TreatRelControl]*_b[WhDisc:Angry]

nlcom [Angry]_b[TreatRelControl]*[WhDisc]_b[Angry] + [Worried]_b[TreatRelControl]*[WhDisc]_b[Worried] + [Disgusted]_b[TreatRelControl]*[WhDisc]_b[Disgusted] + [Afraid]_b[TreatRelControl]*[WhDisc]_b[Afraid]
display .082 / (.082+.081)


*** Christian discrimination

* Race treatment

sem (Disgusted -> CDisc, ) (Angry -> CDisc, ) (Afraid -> CDisc, ) (Worried -> CDisc, ) (PartyID -> CDisc, )  (Ideo -> CDisc, )  (SR -> CDisc, ) (SR -> Disgusted, ) (SR -> Angry, ) (SR -> Afraid, ) (SR -> Worried, ) (BA -> CDisc, ) (BA -> Disgusted, ) (BA -> Angry, ) (BA -> Afraid, ) (BA -> Worried, ) (Attend -> CDisc, ) (Attend -> Disgusted, ) (Attend -> Angry, ) (Attend -> Afraid, ) (Attend -> Worried, ) (Edu -> CDisc, ) (Income -> CDisc, )  (Male -> CDisc, ) (Age -> CDisc, ) (TreatRaceControl -> CDisc, ) (TreatRaceControl -> Disgusted, ) (TreatRaceControl -> Angry, ) (TreatRaceControl -> Afraid, ) (TreatRaceControl -> Worried, ), vce(bootstrap, rep(500) dots(1))  cov( e.Disgusted*e.Angry e.Disgusted*e.Afraid e.Disgusted*e.Worried e.Angry*e.Afraid e.Angry*e.Worried e.Afraid*e.Worried) nocapslatent

estat gof, stats(all)
estat eqgof

nlcom _b[Worried:TreatRaceControl]*_b[CDisc:Worried]
nlcom _b[Disgusted:TreatRaceControl]*_b[CDisc:Disgusted]
nlcom _b[Afraid:TreatRaceControl]*_b[CDisc:Afraid]
nlcom _b[Angry:TreatRaceControl]*_b[CDisc:Angry]

nlcom [Angry]_b[TreatRaceControl]*[CDisc]_b[Angry] + [Worried]_b[TreatRaceControl]*[CDisc]_b[Worried] + [Disgusted]_b[TreatRaceControl]*[CDisc]_b[Disgusted] + [Afraid]_b[TreatRaceControl]*[CDisc]_b[Afraid]
display .01 / (.01 + .017)

* Religion treatment

sem (Disgusted -> CDisc, ) (Angry -> CDisc, ) (Afraid -> CDisc, ) (Worried -> CDisc, ) (PartyID -> CDisc, )  (Ideo -> CDisc, )  (SR -> CDisc, ) (SR -> Disgusted, ) (SR -> Angry, ) (SR -> Afraid, ) (SR -> Worried, ) (BA -> CDisc, ) (BA -> Disgusted, ) (BA -> Angry, ) (BA -> Afraid, ) (BA -> Worried, ) (Attend -> CDisc, ) (Attend -> Disgusted, ) (Attend -> Angry, ) (Attend -> Afraid, ) (Attend -> Worried, ) (Edu -> CDisc, ) (Income -> CDisc, )  (Male -> CDisc, ) (Age -> CDisc, ) (TreatRelControl -> CDisc, ) (TreatRelControl -> Disgusted, ) (TreatRelControl -> Angry, ) (TreatRelControl -> Afraid, ) (TreatRelControl -> Worried, ),  vce(bootstrap, rep(500) dots(1)) cov( e.Disgusted*e.Angry e.Disgusted*e.Afraid e.Disgusted*e.Worried e.Angry*e.Afraid e.Angry*e.Worried e.Afraid*e.Worried) nocapslatent

estat eqgof
estat gof, stats(all)

nlcom _b[Worried:TreatRelControl]*_b[CDisc:Worried]
nlcom _b[Disgusted:TreatRelControl]*_b[CDisc:Disgusted]
nlcom _b[Afraid:TreatRelControl]*_b[CDisc:Afraid]
nlcom _b[Angry:TreatRelControl]*_b[CDisc:Angry]

* total indirect effect
nlcom [Angry]_b[TreatRelControl]*[CDisc]_b[Angry] + [Worried]_b[TreatRelControl]*[CDisc]_b[Worried] + [Disgusted]_b[TreatRelControl]*[CDisc]_b[Disgusted] + [Afraid]_b[TreatRelControl]*[CDisc]_b[Afraid]
display .085 / (.085 + .076)

